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We formulate the basic framework of thermodynamical entropic force cosmology which allows 
variation of the gravitational constant G and the speed of light c. Three different approaches to the 
formulation of the field equations are presented. Some cosmological solutions for each framework are 
given and one of them is tested against combined observational data (supernovae, BAO, and CMB). 
From the fit of the data it is found that the Hawking temperature numerical coefficient 7 is two to 
four orders of magnitude less than usually assumed on the geometrical ground theoretical value of 
0(1) and that it is also compatible with zero. Besides, in the entropic scenario we observationally 
test that the fit of the data is allowed for the speed of light c growing and the gravitational constant 
G diminishing during the evolution of the universe. We also obtain a bound on the variation of c 
to be Ac/c oc 10~® > 0 which is at least one order of magnitude weaker than the quasar spectra 
observational bound. 
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I. INTRODUCTION 

General Relativity is an established theory which ex¬ 
plains the evolution of the universe on large scales [1]. 
Although it is not complete because it contains singular¬ 
ities, it explains the dynamics of the universe in a con¬ 
sistent way. Besides, the current phase of accelerated 
evolution of the universe has been discovered O [3] . In 
order to obtain this accelerated expansion, one has to 
put an extra term, the cosmological constant A or dark 
energy into the Einstein-Friedmann equations. Resulted 
AGDM models [IHZ] are consistent models to explain this 
accelerated expansion, but the observational value of A is 
over 120 orders of magnitude smaller than the value cal¬ 
culated in quantum held theory, where it is interpreted 
as vacuum energy. This motivates cosmologists to look 
for alternative models which can explain the effect BM- 

The relation between the Einstein’s gravity and ther¬ 
modynamics is a puzzle. In the seventies of the twen¬ 
tieth century, Bekenstein and Hawking [ini E] derived 
the laws of black hole thermodynamics which emerged 
to have similar properties as in standard thermodynam¬ 
ics. Jacobson [T2] derived Einstein held equations from 
the hrst law of thermodynamics by assuming the propor¬ 
tionality of the entropy and the horizon area. A more 
extensive work in this direction was made by Verlinde 
and Padmanabhan in Refs. [TSHU]. Verlinde derived 
gravity as an entropic force, which originated in a sys¬ 
tem by the statistical tendency to increase its entropy. 
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He assumed the holographic principle which stated 
that the microscopic degrees of freedom can be repre¬ 
sented holographically on the horizons, and this piece of 
information (or degrees of freedom) can be measured in 
terms of entropy. The approach got criticized on the base 
of neutron experiments though m- 

Recently, the entropic cosmology based on the notion 
of the entropic force was developed in series of papers 
[19H23j and especially it was compared with supernovae 
data in Ref. |24j . However, supernovae test is not very 
strong and so the Ref. [23] got criticized on the base of 
galaxy formation problem (e.g [211 IH|). Basically, the 
idea of entropic cosmology is to add extra entropic force 
terms into the Friedmann equation and the acceleration 
equation. This force is supposed to be responsible for 
the current acceleration as well as for an early exponen¬ 
tial expansion of the universe. It is pertinent to mention 
that the entropic cosmology suggested in these references 
assumes that gravity is still a fundamental force and that 
it includes extra driving force terms or boundary terms in 
the Einstein held equations. This is unlike Verlinde m, 
who considers gravity as an entropic force, but not as a 
fundamental force (see also Refs. [25U39j ). All frame¬ 
works were discussed in detail in Ref. [40] by Visser. 
Entropic cosmology is also related to dynamical vacuum 
energy models which have been discussed and confronted 
with data in Refs. [311445] . 

In this paper we expand entropic cosmology suggested 
in Refs. [TnH23] for the theories with varying physical 
constants: the gravitational constant G and the speed 
of light c. Although [23] is problematic in the context 
galaxy formation test, we use it as a starting point for 
further discussion. We discuss possible consequences of 
such variability onto the entropic force terms and the 
boundary terms. As it has been known for the last 
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fifteen years, varying constants cosmology gSlIlTl was 
proposed as an alternative to inflationary cosmology, be¬ 
cause it can to solve all the cosmological problems (hori¬ 
zon, flatness, and monopole). In the paper we try three 
different approaches to formulate the entropic cosmology 
with varying constants. In section |ll] we present a con¬ 
sistent set of the field equations which describes vary¬ 
ing constants entropic cosmology with general entropic 
force terms. In Section [ITT] we derive the continuity equa¬ 
tion from the first law of thermodynamics and fit general 
entropic terms to the field equations derived in Section 
|TT] using explitic definitions of Bekenstein entropy and 
Hawking temperature. We also discuss the constraints on 
the models which come from the second law of thermody¬ 
namics. In section |TV] we study single-fluid accelerating 
cosmological solutions to the field equations derived in 
Section [Till In Section]^ we derive the entropic force for 
varying constants, define appropriate entropic pressure, 
and modify the continuity and acceleration equations. 
We also determine Friedmann equation and give single¬ 
fluid accelerating cosmological solutions. In Section [Vl| 
we derive gravitational Einstein field equations using the 
heat flow through the horizon to which Bekenstein en¬ 
tropy and the Hawking temperature is assigned. Sections 
VII and VHI| are devoted to observationally testing the 


many-fluid entropic force models with varying constants. 
For this sake the data from supernovae, Baryon Acoustic 
Oscillations (BAO), and Cosmic Microwave Background 
(CMB) is used. In Section IX we give our conclusions. 


II. ENTROPIC FORCE FIELD EQUATIONS 
AND VARYING CONSTANTS 


which will further be used in our paper to various ther¬ 
modynamical scenarios of the evolution of the universe. 
It is clear that (H.3) has dissipative terms in full anal¬ 


ogy to bulk viscosity models. However, if the functions 
/(<) and g{t) are equal and have the value of the A-term 
modified by varying speed of light c(t) i.e. 


fit) = git) = 


Ac^(t) 


(11.4) 


then they give modified varying c and G Einstein field 
equations with the continuity equation as |48j 


0 + 3 - 61 - 1 - 


c^t) 


G(t) (3fc — Aa^) 
■ G(t) ^ 47rG(f)a' 


2 c(t)c(f). 


(H.5) 

which finally reduce to the standard A-CDM equations 
for c and G constant. Another point is the f{t) and g{t) 
terms can also be considered as time-dependent (dynam¬ 
ical) vacuum energy 


III. GRAVITATIONAL THERMODYNAMICS 
AND VARYING CONSTANTS 

In this section we start with basic thermodynamics in 
order to get entropic force varying constants field equa¬ 
tions. We remind that the first law of thermodynamics 
has been widely used to interlink different gravity theo¬ 
ries with thermodynamics P5H551 . Defining the tem¬ 
perature and entropy on the cosmological horizons, one 
can use this law of thermodynamics for the whole uni¬ 
verse 


The main idea of our consideration is to follow Refs, 
[m HOI m] (assuming homogeneous Friedmann geome¬ 
try) and generalize field equations which contain the en¬ 
tropic force terms f{t) and g{t) onto the case of varying 
speed of light c and varying Newton gravitational con¬ 
stant G theories. It is easy to realize that the modified 
Einstein equations can be written down as follows 



d 
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87rG(f) kc^(t) 

-^^p-—^ + fit): 


47rG(f) 


3p 

+ (t) 


+ git)- 


(ii.i) 

(H.2) 


dE + pdV = TdS, (HI.l) 

where dA, dV, and dS describe changes in the internal 
energy E, the volume V, and the entropy S, while T is 
the temperature, and p is the pressure. The volume of 
the universe contained in a sphere of the proper radius 
r* = a(t)r (r is the comoving radius and a{t) is the scale 
factor) is 


V{t) 


4 

3 


3 3 

Tra r 


(HI.2) 


We have 


V{t) = = W{t)H{t) (HI.3) 


In fact, the functions f{t) and g(t) in general play the 
role analogous to bulk viscosity (cf. Refs. [50-69] o f the 
paper by Komatsu et. al. [19) 1 and this is why from dilTl ) 
and (H.21 one obtains the modified continuity equation 


where dot represents the derivative with respect to time 
and the Hubble function is H{t) = d/a. The internal 
energy E and the energy density e(t) of the universe are 
related by 


P + 3H 



_L _ o kc{t)c{t) 

G{t) 47rG(f)a2(f) 


3H 

47rG(f) 


git) - fit) 


m 

2H 


(H.3) 


E{t) = e{t)V{t), eft) = p{t)c^{t), (HI.4) 

where p is the mass density of the universe. 

Now we generalize the Hawking temperature T m 
and Bekenstein entropy S HO] of the (time-dependent) 
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Hubble horizon at r = = r/j (t) onto the varying c and 

G theories as follows 


T = 


'yhc{t) 


2-KkBrh{t) ’ 


Q 

4fi 


c\t)A(t) 

G{t) 


(IIL5) 

(IIL 6 ) 


Here A{t) = 47 rr^(t) is the horizon area, h is the Planck 
constant, ks is the Boltzmann constant, and 7 is an ar¬ 
bitrary, dimensionless, and non-negative theoretical pa¬ 
rameter of the order of unity 0 ( 1 ) which is usually taken 
to be ^ or i (THl HOI 111]• In fact, 7 can be related 
to a corresponding screen or boundary of the universe 
to define the temperature and the entropy on that pre¬ 
ferred screen. Here the screen will be the Hubble horizon 


i.e. the sphere of the radius r^,. Dividing (HI.l) by time 
differential dt, we have 


dE dV dS 

dt ^ dt dt' 


(HL7) 


which after applying (HL2) and (HL4) gives 

P 


E+pV = 


„c(t) , 

P + 2—ttP + 3- P + 
c{t) r 




Vc^it). 


From (HL5) and (HL 6 ) we have 


TS = 

2G{t) ^ 


^ 

c{t)^ rh G{t) 


(HL 8 ) 


(HL9) 


By using (HI.71, (HI. 8 ), and (HI.9) we get the modified 


continuity equation as follows 

8TTG{t) [ c{t) G{t) 


(HI.IO) 


- 2 ^ 

H 


which for A: = 0 case reduces to 


r_h 

Th 


c{t) 

c{t) 


H 

H' 


(HI. 14) 


In this section we restrict ourselves to fc = 0 case in order 
to get the general functions f{t) and g{t). 

In order to constrain possible sets of varying constant 
models we can apply the second law of thermodynam¬ 
ics according to which the entropy of the universe re¬ 
mains constant (adiabatic expansion) or increases (non- 
adiabatic expansion) 


dS 

dt 


> 0 . 


(HI.15) 


In fact, (HI.91 gives the condition 


p Cji) _ G^ _ ^ ^ 
C{t) G{t) - Th 



^ . (HI.I 6 ) 


or 


5'S(-((M>2*=2(“-2) (IIU7) 

c{t) G{t) El \d a) 

which for c = G = 0 just says that the Hubble horizon 
must increa se Vh > 0. For G[t) = 0, and by using (HI. 11) 
and (HI. 16), we have 


c {t) > biH 5 , 
and for c{t) = 0, we have 

G{t) < 

where bi and 62 are constants. 


(HI.18) 

(HI.19) 


IV. GRAVITATIONAL THERMODYNAMICS - 
COSMOLOGICAL SOLUTIONS 


where we have used the explicit definition of the Hubble 
horizon modified to varying speed of light models [24] 


rh{t) = 


c(A) 

H{ty 


(HI.ll) 


If we introduced the non-zero spatial curvature fc = ±1, 
then we would have to apply the entropy and the tem¬ 
perature of the apparent horizon which reads 


TA = 


it) 


V" + aHt) 


(HI.12) 


Simple calculations give that 


Using the generalized continuity equation (HI.IOI one 


is able to fit the functions f{t) and g{t) from a gen¬ 
eral varying constants entropic force continuity equation 
(H.3| as follows 


fit) = 


g{t) = 777 ^ + 1(51^- 


.c{t) G{t) 


c{t) G{t) ^ 

47rG(t) /G(t) _ „cii)\ 
3H [ G{t) c{t) i 


H 


(IV.l) 


(IV.2) 


Having given f{t) and git) one is able to write down the 
equations (H.l[) and (H.2) as follows 


rA 


Hr\ 


kA 


77-— +- l--rl 

7 c V a 


(HI.13) 


+ (IV.3) 
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A. G varying models only: G{i) = GoO,'^', 
q, Go = const., c{t) — 0. 
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- = - 
a 


47rG(t) 


77 - 2 ^ c{t) 
z{t) 


H 


3p 

l-27^^ G(t) 

W) 


H, (IV.4) 


which form a consistent set together with Eq. (III.10). 


While fitting the functions f{t) and g{t) we set fc = 0. 
If we were to investigate fc = ±1 models then the the 
temperature T (III.5) and the entro py S ( III.6) should 
be defined on the apparent horizon (III. 12). A different 


choice of /(t) and g(t) which is consistent with (III.10) 
would be for example as follows 


m = 0 , 


+ 


47rG(t) 

ZH 


Gjt) 

G{t) 


- 2 


c{t) 

c{t) 


P- 


(1V.5) 


(IV.6) 


Howev er, both choices (IV.l )-( |rW2 ) and (IV.5) and 


(IV.6) do not allow for a constant term like the cosmo¬ 
logical constant (unless one fine-tunes H = const.) so 
that an alternative choice which fulfills this requirement 
would be 


/(t) = 7ij" -f Ai, 
g(t) = -fH'^ + Ki + 


(IV. 7) 


7 ( .cW _ ^ 
2 I c{t) G{t) 


H 


47rG(t) / G(t) _ 2 ^ ] . 
3H I G{t) c{t) ] 


(IV.8) 


where Ki is a constant acting on the same footing as 
the cosmological A-term in standard A-CDM cosmology 
securing the model with respect to structure formation 
tests (cf. [HIIMIS]). 

There is a full analogy of varying constants generalized 
equations (IV.3), (IV.4), and (III. 10) with the entropic 


force equation given in |24) when one applies the specific 
ansatze for varying c and G: 


c{t) = CoO" 

with n,q = const, 
or G{t)/G{t) = qH. 


and G{t) = Go 


(IV.9) 


which gives c{t)/c{t) = nH 

It is worth emphasizing that 
our ansatze should be c(t) = Co{a/ao)'^ and G{t) = 
Go{a/ao)‘‘ but the standard approach nowadays 
picks up oq = 1 [S3]. _ 


We note that the application of the ansatze (1V.9) to 


the growing entropy requirement (111.17) gives the bound 
that 


n > —2/5 


and 


q<2. 


(IV.IO) 


As we shall see in Section Vlll (or Table 1) these limits 
are in agreement with the observational values we have 
obtained. They also allow the Newtonian limit c —^ 00 
(n —>■ 00 ) or G —0 (q —>■ — 00 ) [S 0 | . 

The cosmological solutions of the set of varying con¬ 
stants Eqs. (IV.3), (IV.4) and (III. 10) are given below. 


Defining the barotropic index equation of state param¬ 
eter w by using the barotropic equation of state, p = wpc^ 
for varying G = GoOp, we can integrate the continuity 
equation ( 111 . 10 ) to get 


P = 


{ T ]{1 


(IV. 11) 


where po is a constant with the dimension of mass den¬ 
sity, Go the gravitational constant, and Hq the Hubble 
param eter. Using ( |IV.3 ) and ( IV.4) and then multiplying 
(IV.3) by (I -I- 3w) and (IV.4) by 2 we get 


H = 


+ (IV.12) 


or (using the fact that GjG = qH, d/a = H + H"^) one 
has 


H = -wH^, 


where 


w = i [3{w -h 1) (1 - 7 ) - (1 - 2-f)q\. 


(IV.13) 


(IV. 14) 


In a [in] i.e. 


The Eq. (IV.13) solves easily using a new variable N = 

(1V.15) 


dH \ dH (dH \ dt _ 
-r- o, = —— = —— —a = —wH, 

da I dN \ dt I da 


which integrates to give 

H = Hoa- 


(1V.16) 


where Hq is constant. The solution of (IV.16) is 
a{t) = [Ho{t-to)]^ , 


(IV. 17) 


where to is constant. Bearing in mind the value 
of (IV. 14), we can easily conclude that without en¬ 


tropic terms the solution (IV. 17) corresponds to a stan¬ 


dard barotropic fluid Friedmann evolution a(t) cx (t — 
io)(2/3(i«-i-i))^ The scale factor for radiation, matter and 
vacuum (cosmological constant) dominated eras reads as 


a{t) oc 


[Ho{t — to)] ( 4 -‘!)+ 2 v( 9 - 2 ) (radiation) 


[Ho{t — to)] (3-9) + ( 29 - 3)7 ; y; = 0, (duSt) 

[i7o(t — to)] ( 27 - 1)9 = — 1 . (vacuum) 

(1V.18) 

The solution ( IV.18[ ) shows that in varying G entropic 
cosmology even dust (w = 0 ) can drive acceleration of 
the universe provided 

(3-q) + (29-3)7<2 . 


(1V.19) 
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On the other hand, the solution which includes A—term 
{w = —1) drives acceleration for (27 — 1) q < 2. There is 
an interesting check of these formulas for the case when 
one takes the Hawking temperature parameter 7 = 1, in 
all three cases (radiation, dust, vacuum) the conditions 
for accelerated expansion fall into one relation q < 2. In 
fact, this limit is very special which can be seen from Eq. 

) in which the terms involving cancel and lead 
to empty universe (g = 0 ) so that it is no wonder that 
the acceleration does not depend on the barotropic index 
parameter w. Finally, we conclude that in all these cases 
the entropic terms and the varying constants can play 
the role of dark energy. 

One may also consider a more than one component 
models i.e. the models which allow matter, radiation as 
well as other cosmological fluids of negative pressure like 
the cosmological constant which give a turning point of 
the evolution compatible with current observational data 
(early-time deceleration and late-time acceleration). We 
will consider such models numerically in Section |VIII 
where we test these models with observational data. 


(IV. 


B. c varying models only: c{t) = Cofl"; co,n = const., 
G{t) = 0 


The solution of the continuity equation (III. 10 1 for 
varying c is 




c{t) 

Co 




- 27 ' 


, (1V.20) 


where again po is a constant with the dimension of mass 
density, Co the velo city o f light , and Hq the Hubble pa¬ 
rameter. Applying (IV.3|) and (IV.4), we have 


Z c 2 


or 


H = -wH\ 


where 


^ = 2 + w) (1 - 7 ) - n (77 - 2 )]. 


The solution of (IV.22) is 


H = Hoa- 


(IV.21) 

(IV.22) 

(IV.23) 

(IV.24) 


where Hq is constant. Finally, the solution of ( |IV.24 ) for 
the scale factor gives 


i{t) = u)™ [Ho{t- to)]’ 


(IV.25) 


where to is constant. For radiation, dust and vacuum we 
have, respectively 

[Ho{t — to)] ( 4 + 2 ")-( 4 + 7 r,)^ .y; _ 1 ^ (radiation) 
a(t) oc [i7o(t — to)] (3+2")-(3+7 ")t ; 10 = 0, (dust) 

[i7o(t — to)] ; uu = — 1 . (vacuum) 

(IV.26) 


For these three cases, one derives inflation provided 
(4 -I- 2n) — (4 -I- 7n)7 < 2, (radiation) 

(3 -I- 2n) — (3 -I- 7n)-j < 2, (dust) 

(2 — 77)71 < 2, (vacuum) 

and the entropic force terms play the role of dark energy 
which can be responsible for the current acceleration of 
the universe. As in the previous subsection, here also af¬ 
ter taking the Hawking temperature parameter 7 = 1, in 
all three cases (radiation, dust, vacuum) the conditions 
for accelerated expansion fall into one relation n > —2/5, 
but this is also a special empty universe limit of Eq. 

( |!v3t . 

As in the previous subsection one may also consider 
a more than one component models - the matter we 
deal with numerically in Section |VHI| We would like 
to emphasize again that here we have presented one- 
component solutions only while in Section [VH| we will be 
studying multi-component models which allow the tran¬ 
sition from deceleration to acceleration. 


V. ENTROPIC PRESSURE MODIFIED 
EQUATIONS 


In this section we start with the formal definition of 
the entropic force as given in [M HOI mi- We a ssume 
that the temperature and entropy are given by (HI.5) 


and (HI .61 and use the definition of the entropic force 


F = -T^. 
dr 


(V.l) 


We calculate the entropic force on the horizon r = rh{t) 
by taking 

dS/drt = S/r-h (V. 2 ) 


to obtain 


2G{t) 


cyd) _ G{t) _ 

c(t) _ H 
clt) H 


(V.3) 


For c = G = 0 this formula reduces to the value obtained 
in Ref. [19]: F = 'y{c^/G) which is presumably the value 
of maximum tension in general relativity |54H56j . It has 
been shown in Ref. m that ( |V.3| ) may recover infinite 
tension thus violating the so-called Maximum Tension 
Principle [54] in the framework of varying constants the¬ 
ories. 

Now, we define the entropic pressure pe, as the en¬ 
tropic force per unit area A, and use (HI.ll) to get 


7C^(t)i7^ 
^ “ 87rG(t) 


cyCO G{t) r >H 

^cjt) G(t) 

c(t) _ H 
c(t) H 


(V.4) 
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Out of the set of initial equations (|II.l 
of them are independent. On the ot 


ner 


-( |IL3 ) onl y two 
rand, only (II.2) 


(acceleration equation) and (11.3) (continuity e quat ion) 
contain the pressure. This is why while having (V.4) we 
will define the effective pressure 


Peff =P+PE 


(V.5) 


and then write down the continuity equation (II.3) as 

(V.6) 


p + 3i/ I p + 


or 


p + iH + 

87rG(t) 



/ G{t) 

p 

\ , m 

cHt), 

I G{t) 

»c(i) 

Git) 

Git) 


p = 0, 


(V.7) 


^ H 


c(t) _ H 
c(t) H 


and the acceleration equation (II.2) as 
47rG(t) 


a 

a 


^Peff 

c^t) 


(V.8) 


or 


47 rG(t) f ip 

P~^ „2 


it) 




_ G{t) 
°c(t) G(t) 


- 2^ 


c{t) _ H 

c{t) H 


(V.9) 

In order to solve the continuity equation ( V.8[) we have to 


put f{t) = 0 in the Friedmann equation. Alternatively, 
we see by comparing (V.8) and (II.3) for fc = 0, that we 


need to put /(t) = 0. We then obtain the simplest form 
of the Friedmann equation to use 


87rG(t) 


(V.IO) 


By using ( |V.9 ) and (V.IO), we get for varying c{t) = 
Co a” and G(t) = Goa'*: 


where 


Bi^) + 


Bi = —3(1 + w) +2j 
B 2 = 3(1 + re) — 57 


(V.ll) 

(V.12) 

(V.13) 


The cosmological solutions are obtained below. We 
consider two cases. 


A. G varying models only: G{t) V 0 and c(t) = 0; 
g V 0, n = 0. 


The Eq. ( |V.11[ ) reduces to 
2 


(-) 

or we can write 


-if = 


(?) 


or 


where 




H = -WH^, 


(V.14) 


(V.15) 


(V.16) 




Solving (V.16) for the Hubble parameter, we have 
H = Hoa-^, 


(V.18) 


where Hq is a constant of integration. Solving (V.18) for 
the scale factor a{t), one gets 


a(t) = W^ [Hoit-to)]'^ . 


(V.19) 


B. c varying models only: G{t) — 0 and c{t) V 0; 
g = 0, n V 0 


From (V.ll) we obtain 



B-y T 2ti 





(V.20) 


Following the same procedure as in the subsection A, one 
can find the Hubble parameter and the scale factor for 
varying c as: 


H = Hoa-^ (V.21) 


and 

a{t)=xi[Ho{t-to)]^ , (V.22) 

where, Hq and to are real constants and X is given by 


nB 2 (Hi+2n)^\ i?i+2n\ 

_ + j + _ j . 

(V.23) 

Both of the above cases have the same non-varying con¬ 
stants limit (n —>■ 0 or G —>■ 0) of W = By 12. 


















































7 


VI. GRAVITATIONAL THERMODYNAMICS - 
HORIZON HEAT FLOW 


or after explicitly using (111.12) we get a generalized ac¬ 
celeration equation 


In this section we use yet another approach to de¬ 
rive entropic cosmology which is based on the applica¬ 
tion of the idea that one can get gravitational Einstein 
field equations using the heat flow through the horizon to 
which Bekenstein entropy (III.6) and Hawking tempera¬ 
ture (III.51 (with 7=1) are assigned. 

The heat flow dQ out through the horizon is given by 
the change of energy dE inside the apparent horizon and 
relates to the flow of entropy TdS as follows 


dQ = TdS = -dE. 


(VI.1) 


If the matter inside the horizon has the form of a perfect 
fluid and c is not varying, then the heat flow through the 
horizon over the period of time dt is m 


dQ rr.dS P 2 / , Pn 

-77=^37= = 47rr^(p+ 


(VI.2) 


However, in our case c is varying in time and we have to 
take this into account while calculating the flow so that 
bearing in mind that the mass element is dM we have 
the energy through the horizon as 


-dE = c^dM + 2Mcdc + pdV. (VI.3) 


The mass element flow is 


dM = A{vdt)g = dV g, (VI.4) 

where vdt = s is the distance travelled by the fluid ele¬ 
ment, V is the velocity of the volume element, and dV is 
the volume element. The velocity of a fluid element can 
be related to the Hubble law of expansion 


V = EIta 


so that (VI.41 can be written down as 


dM = AEIrAgdt. 


(VI.5) 

(VI.6) 


We assume that the speed of light is the function of the 
volume through the scale factor i.e. c = c(V) and since 
a oc then c = c(a) |SH]. We have 


dc 1 1 dc 

dV ^ 3 y2/3 da 


(VI. 7) 


and besides by putting M = Vgin (VL3) we get 


.rn 2 .Tr / P 2 a dc\ 
-dE = c^dV { g+-^ + ■ 

C o C CLd j 


(VI.8) 


Using ( VI.6), (VI.8 ) and (HL9) (replacing ry, by r^) one 
has from (VI.l) 


AtttaH 



2 a dc\ 
3^c da) 


2G I c ^ G 


(VI.9) 


i7 = -4.G(p+J)-bM^-gU 


8ttG c p 1 k(? 


G 


3 cH ■ 


which for c = G = 0 gives the Eq. (A6) of Ref. [24] . 

In order to get the Friedmann equation, we have to 
use the continuity equation (VIL22) but for adiabatic 
expansion {dS = 0) to obtain 


H 




p 2 c 


P- 


By using Eq. (VI. 10) in ( jVI.ll ), we have 




(VI.ll) 


1 kc^ It G 

2lM I c “ G 


2H 


(VI.12) 


After integrating (VI.12) one obtains a generalized Fried¬ 
mann equation 


= ^ f G{t)pdt + 




kc^ 




(VI.13) 


For k = 0 [ta Th = c/ El) by taking the ansatz of the 
form 


c{t)=co[H{t)r, (VL14) 

Co = const., m = const, (or c{t) = Co{H/H q)"^, Hq = 
const.; similar ansatz c{t) = d{t) was used in Ref. [59]), 
for varying c only (i.e. for G = 0), we have the following 
equations 


o SttG 5m 
^ = -^P + 


(VI.15) 


d AttG , 3p, 5m ,3m bnnS. ■ H 

(VI.16) 


p + S-ff I p + 


\t) 


+ 2m—p = 0, 


(VI. 17) 


where solely K is the constant of integration which can 
be interpreted as the cosmological constant Ac^/3 [55] 
provided m = 0 (cf. the discussion in Section [IT] and the 



































formula (II.5)). For small m one may expand c{t) given 
by (VI. 14) in Taylor series 

(VI.18) 


c{t) = Co [Hit)]"^ 

1 + m\nH{t) 


= Co 




and use numerical procedures to calculate the conse¬ 
quences of variability of the speed of light but we keep 
this beyond the scope of the paper. 

Our set of equations ( VI.15| )-( VI.17) contains two ef¬ 
fects: the entropic force contribution K as well as many 
new terms related to variability of c (all the terms which 
involve the parameter m). In fact, there are as many as 
four such latter terms in the equation (VI. 16) (including 
a cross-term with K) and each of them may play the role 
in accelerating the universe instead of K-ieim.. 

For K = 0 case one can easily solve for the Hubble 
parameter H and the scale factor a for varying c models 
as 


and 


3(l + »^) 

H = Hoa 2(1+™), 


2(1 + ™) 

a(t) (x [Ho{t - to)] 


(VI.19) 


(VI.20) 


where Hq is the constant of integration. Besides, the 
continuity equation solves by 


3(1 + ™) 

p = poa i+™ , 


(VI.21) 


where po is a constant with the dimension of mass density. 
The solutions for K ^ 0 can be found numerically, but 
we do not present them here. 


VII. OBSERVATIONAL PARAMETERS 

In this section we will try to give some more quanti¬ 
tative information about our approach, by applying our 
model to observational data. We will leave the single¬ 
fluid approach we have considered in past section, to 
move to the more realistic case of a multi-fluid scenario. 
We will take into account the components which make 
up the total mass density p, i.e., radiation pr, mat¬ 
ter Pm, and some unknown vacuum energy component 
Py (which can also be the cosmological constant). We 
take the model with f{t) and g{t) given by (IV.l) and 


(IV.2) which do not contain the constant Ai-term as 
in ( IV.7| ) and (IV.8). However, we will get this con¬ 
stant term effectively as the energy density of vacuum 
PA = {Ac^)/{8ttG). Wit h these assumptions we can write 
the continuity equation (HI.10) as 




E- 


c^t) 


= (VH.22) 


7 


+ T^E^* 

1 — 7 


At) cm H 
, c{t) G{t) H 


where summation on i runs for radiation, matter and 


dark energy. From (VH.22) one can easily check that one 


can separate the contribution of the three fluids obtaining 
three separate continuity equations; if we introduce the 
equation of state parameter Wi from pi = Wi pi, then the 
continuity equations have the form: 


Pi + iHpi (1 -I- Wi) = -2-pi 
c 


(VH.23) 


1-7 


Pi 


.c{t) G{t) 


c(t) G{t) 



where, of course, Wi = 0 for matter, Wi = 1/3 for radia¬ 
tion, and Wi = —1 for vacuum. From them, it is easy to 
check that, on one hand, no interaction term is present, 
in the way of exchanging energy among the fluids; but, 
on the other, entropic forces and the varying constants 
influence the behaviour of the fluids, by the same amount 
and separately. Thus, for each of them, a separate con¬ 
tinuity equation holds, and we never have any violation 
of the mass-energy conservation law. 


The solution for each fluid from (VH.23) can be easily 
found; once we use our ansatze, c = cga" and G = Goo'^, 
we have: 


Pi = 


Po 


2 + 

HA 


iJTd/a/f(7.n,9) 


(VH.24) 


where, as usual, H is the Hubble function. Ho the Hub¬ 
ble constant, a the scale factor, and fA , n, q) are general 
functions obtained by solving (VH.23). When consider¬ 


ing only a varying c, these functions are: 


fii'1,1'1) = -3 


l+Wi + 


n(2 - 77 ) 

3(1 - 7 ) J 


while for a varying G we have: 


/f(7,g) = -3 


1 -h Wi -\- 




3(1-7) 


(VH.25) 


(VH.26) 


We can note down that the main changes to the equa¬ 
tion of state parameters come for the varying constant 
assumptions: in the limit of n —>■ 0 and (? —>■ 0, we re¬ 


cover the usual behaviours, a~'^ for matter, for radi¬ 
ation, and constancy for the vacuum. But still we have 
some dynamical effects on the densities from the entropic 

27 

forces, through the H ^-1 term. Thus, even in the case 
of no-varying constant, the entropic forces make the vac¬ 
uum dynamical. 


Starting from the Friedmann equation (IV.3), after 


some simple algebra, we can write the Hubble function 
H, which we explicitly need for observational fitting. In 
the case of varying c it will be 


E^=(^ 


\Ho 


a 


2,0 

— - nJz 


“1-7 


{l,n) 


T-1 

27-1 


(VH.27) 
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while for varying G it will be 


= 


H 


^ 1 — 7 


( 7 , 9 ) 


7-1 

27-1 


(VIL28) 


We have defined the dimensionless density parameters as 


= 


SnGoPifi 

-JhT' 


(VIL29) 


where, Gq is the current value of Newton’s gravitational 
constant. Finally, in order to check if our model (VIL22) 


allows a transition from deceleration to acceleration dur¬ 
ing the evolution of the universe at some redshift z in a 
similar way to a “pure” ACDM model, we have looked 
at the deceleration parameter, defined as: 


with H{z) = HqE{z) (following [50], we assume i7o = 70 
km/s Mpc“^), Co the speed of light here and now, and 0^ 
the vector of cosmological parameters. The total vector 
6 will include 6c and the other fitting parameters, which 
in this case are: a and /3, which characterize the stretch- 
luminosity and color-luminosity relationships; and the 
nuisance parameter Mb, expressed as a step function 
of two more parameters, and A^: 


Mb = 


M], 

M\, 


< 


Q, 


if Alg^clla 
Am otherwise. 

(VIII.34) 

Further details about this choice are given in Ref. m- 
The formula (VIII.35) stands for th e constant c cases; 
when c is varying according to (IV.9), it is modified into 

Emm 


9(^) 


(l + z) dH^{z) 
2H'^{z) dz 


(VII.30) 


Dl{z,9c) 


cq 


(1 + ^) 



{1 + z')-^ 
E(z',9c) 


dz' . 


(VIII.35) 


where the cosmological redshift is given by l-|-z = 1/a. 


B. Baryon Acoustic Oscillations 


VIII. DATA ANALYSIS 

The analysis has involved the largest updated set of 
cosmological data available so far, and includes: Type 
la Supernovae (SNela); Baryon Acoustic Oscillations 
(BAO); Cosmic Microwave Background (CMB); and a 
prior on the Hubble constant parameter, E[q. 

A. Type la Supernovae 


The X%AO Baryon Acoustic Oscillations (BAO) is 
defined as 

xIao = ■ CsAo • > (VIII.36) 

where the quantity can be different depending on 

the considered survey. We used data from the Wig- 
gleZ Dark Energy Survey m, evaluated at redshifts 
z = {0.44,0.6,0.73}, and given in Table 1 of [52]; 
this case the quantities to be considered are the acoustic 
parameter 


We used the SNela (Supernovae Type la) data from 
the JLA (Joint-Light-curve Analysis) compilation [60] . 
This set is made of 740 SNela obtained by the SDSS-II 
(Sloan Digital Sky Survey) and SNLS (Supenovae Legacy 
Survey) collaboration, covering a redshift range 0.01 < 
z < 1.39. The Xsn defined as 

xIn = • Csl, ■ Ajzs^ , (VIII.31) 


A(z, 9c) = 100VO,„ ^^ Dv{z,6c) ^ 


CqZ 


37) 


and the Alcock-Paczynski distortion parameter 
Da{z,9c)H{z,6c) 


E{z,9c) = {1 + z)- 


co 


(VIII.38) 


where. Da is the angular diameter distance 


with Ajz^^ = , the difference between 

the observed and the theoretical value of the observable 
quantity ; and Csn the total covariance matrix (for 
a discussion about all the terms involved in its deriva¬ 
tion, see m)- I^or JLA, the observed quantity will be 
the predicted distance modulus of the SNela, p,, given the 
cosmological model and two other quantities, the stretch 
(a measure of the shape of the SNela light-curve) and the 
color. It will read 

p{9) = 5\og,o[DL{z,9c)]- aXi+ pc + Mb , (VIIL32) 
where Db is the luminosity distance 


Da{z,9c) 


Co 1 dz' 

Ho 1+^ Jo E{z',ec) ’ 


(VIIL39) 


and Dv is a combination of the physical angular-diameter 
distance Da (tangential separation) and Hubble param¬ 
eter H(z) (radial separation) defined as 


Dv{z,9c) = 


When 


Eqs. ( IVHL37D 
be changed into 


dealing 
(|VIH.38 

EMSfr 


{l + z)^D\{z,9c) 
with 


CqZ 


-I 1/3 


H{z,9c)\ 

(VHL40) 
varying c, 


38) (VHL39) (VHL40) have to 


Dl{z,9c) 


1^(1+ .) 

Ho 


r dz' 

Jo E{z',ec) ’ 


(VHL33) 


A(z, 9c) = lOOV^mh^ ’ (VHL4I) 

Co(l-l-z) "z 
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z a 

FIG. 1. Transition from deceleration to acceleration for the model ( |VII.22[ ). Solid line is for standard ACDM model; dashed 
line for varying-G-entropic-ACDM model; dotted line is for varying-c-entropic-ACDM model. The values of the parameters 
are taken from Table I and from the Planck data (ACDM). 


F(z, 0c) = (1 + z) 


DA{z,e,)H{z,e,) 

Co(l + z)-" 


(VIII.42) 


Da{z,6c) 


cq 1 
Hol + z 



(1 + z')"" 
E{z',e^) 


dz' 


(VIII.43) 


Dv{z, 9c) = 


(1 + z)‘^D\{z, 9c) 


co(l 




-I 1/3 


H{z,9c) 


(VIII.44) 

We have also considered the data from SDSS-III Baryon 
Oscillation Spectroscopic Survey (BOSS) DRlO-11, de¬ 
scribed in [63l [64] • Data are expressed as 


Dv{z = 0.32) 


(1264 ± 25) 

ri (Zd) 


(VIII.45) 


and 


Dy(z = 0.57) = (2056 ±20)^4^, (VIII.46) 

s \^d) 


where rs{zd) is the sound horizon evaluated at the drag- 
ging redshift Zd, and rE‘^{zd) is the same sound horizon 
but calculated for a given fiducial cosmological model 
used, being equal to 149.28 Mpc [HSIIM]- The redshift of 
the drag epoch is well approximated by [SS] 


Zd 


1291(0^ h2)0-251 

1 -b 0.659(0^ /i2)0.828 


[l + bii^bh^)^^] (VIII.47) 


where 


with the sound speed 

Csiz) = - (VIII.50) 

y3(l -|- Rb (1 -b z)~^) 

and 

Rb = 3150006 {Tcmb/2.7)~'^ , (VIII.51) 


with Tqmb = 2.726 K. 

We have also added data points from Quasar-Lyman 
a Forest from SDSS-III BOSS DRll [5^ : 


Da{z = 2.36) 
rs{zd) 

_Co_ 

H{z = 2M)rs{zd) 


10.8 ±0.4 
9.0 ±0.3. 


(VIII.52) 


When working with varying c models, of course, we will 
have to change Da and Dy as described above, and also 
the sound horizon, through the definition of the sound 
speed, Eq. (VIII.50), which now will be 


Cs{z) 


co(l±z)-" 

y^3(i±i?6(i±2^)~^) 


(VIII.53) 


Thus, we will have three different contributions to Xbao^ 
e-g-. XwzggieZ^xloss^xlyman^ depending on the data 
sets we consider. 


C. Cosmic Microwave Background 


bi = 0.313(fl„ /i2)-o-4i9 [1 + 0.607(D„ , 

62 = 0.238(fl„ /i2)0-223. (VIII.48) 

And the sound horizon is defined as: 


The XcMB foi' Cosmic Microwave Background (CMB) 
is defined as 

xImb = , (VIII.54) 


rs 



Cs(z') 

H{z',9c) 


dz ', 


(VIII.49) 


where is a vector of quantities taken from 

where Planck first data release is analyzed in order to 
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give a set of quantities which efficiently summarize the in¬ 
formation contained in the full power spectrum (at least, 
for the cosmological background), and can thus be used 
in alternative to the latter [55] • The quantities are the 
CMB shift parameters: 

* Co 

and the baryonic density parameter, Again, Vg 

is the comoving sound horizon, but evaluated at the 
photon-decoupling redshift z*, given by the fitting for¬ 
mula [55] : 


Z* = 1048 [1 -h 0.00124(flb/l2)-0.738j ^ , 

(VIIL56) 

with 


2\-0.238 


9i = 


92 = 


0 . 0783 (^ 6 / 12 ) 

1-h 39.5(fl6/l2)-0.763 

0.560 


l-h21.1(fl6/l2)l-81 ’ 

while r is the comoving distance defined as: 

dz' 


(VIII.57) 




When considering varying c models, again, the sound 
horizon will change as described above, and the comoving 
distance will be iniiis2] 


r(z,ec) = ^ 


(1 + z')- 


dz' 


i7o7o E{z',e^) 

and the shift parameter R will become 


(VIII.59) 


R{e,) ^ (VIII.60) 


Moreover, we have added a gaussian prior on the Hub¬ 
ble constant, Hq 


Xho 


{Hq - 69.6)2 

0.072^^ 


(VIII.61) 


derived from HQ]. 

Thus, the total Xrot ^e the sum of: 

XsN^Xwiggiez^xloss^xlyman^XcMB^XH„- We min¬ 
imize Xrot using the Markov Chain Monte Carlo 
(MCMC) method. 

Finally, we should make a few comments about 
the parameters which will be constrained. The 
total parameters vector Of. will be equal to 
{nm,^b,h,q,j,a /3 Mg, Am} when considering the 
varying G cases, and {Xtrn,^b,h,n,'y,a P Mg, Am} 
when considering the varying c ones. The actual obser- 
vationally fitted components of this vector are given in 
Table H 


The parameter h is defined in a standard way by 
Hq = 100/i. The density parameters entering H{z) are 
nm,nr,Xly', assuming zero spatial curvature, we can ex¬ 
press Xly = 1 — 7 — Xlm — Xly, in order to ensure the 
condition E{z = 0) = 1. Moreover, the radiation density 
parameter Xly will be defined HH as the sum of photons 
and relativistic neutrinos 

riy = Q^{l + 0.227lMeff), (VIII.62) 

where = 2.469 x 10-^ h-^ for Tcmb = 2.726 K; 
and the number of relativistic neutrinos is assumed to 
he Meff = 3.046. 


D. Results 


Our main result is presented in Fig. First novelty 
is that we have found the observational bounds on the 
Hawking temperature coefficient 7 which (on the theo¬ 
retical basis) was usually taken to of order of unity 0(1). 
Our evaluation gives that it should be of the order of 
10“2 — 10“^. This difference is not unexpected, because 
the 0(1) estimation was based on purely theoretical con¬ 
siderations, with no previous connection to data. Now, 
we show that observations are not consistent with such 
large values of 7. Instead, it is at least two orders of 
magnitude less. Thus, the entropic force in the model 
we have considered gives only a small contribution. Sim¬ 
ilar results were obtained in [mill]. Another novelty is 
the bound on the variability of the speed of light c and 
the gravitational constant G. According to them, in the 
entropic scenario we have investigated both G (Fig. 
left panel) and c (Fig. right panel) should be increas¬ 
ing with the evolution of the universe. Bearing in mind 
that the speed of light is related to the inverse of the fine 
structure constant defined as 


a = 


he 


(VHI.63) 


where e is the electron charge and h is the Planck con¬ 
stant, by using (IV.91 and (VHI.63) one has 


Ac Aa Aa n 

c a a 10’ 


(VHI.64) 


then one can derive from Table [^ and Fig. [^ that the 
change in c and so in a (Aa/ao) from our fit is ~ 10“® 
in a redshift range [1; 2] (n = 4.9 • 10“^ > 0), while other 
observational bounds, in the same range (see table H of 
Ref. [72] which is based on [75H75] 1. give Aa/ao ~ I0“®. 
But still our estimation is compatible with other cosmo¬ 
logical constraints, as the ones derived from CMB Planck 
first release, see m- Moreover, recent observations show 
that both positive and negative values of n are possible 
(the so-called a—dipole [77]). 

Finally, we can enumerate some general conclusions as 
follows: 
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• the entropic scenario plus varying c and/or G is 
quite indistinguishable from a pure-ACDM model, 
that is why we call it an entropic-ACDM model. 
Present data is still unable to differentiate between 
the two scenarios; 

• the model obtained {entropic-ACDM cosmology) is 
a variation of the exchange of energy between vac¬ 
uum and matter model studied in Refs. gmum. 

• the best fit for the value of the Hawking temper¬ 
ature coefficient 7 is quite different from the the¬ 
oretical values used in literature, i.e. 7 = 3/(27r) 
or 1/2; it should be pointed out that other consid¬ 
ered entropic scenarios have the values of 0(1) (e.g. 

my, 

• the model with small values of the parameter 7 is 
equivalent to a dynamical vacuum model with small 
variation of the vacuum energy studied in Refs. [411 

mi; 

• the value for 7 is compatible with zero since we were 
able to put only an upper limit to it. This would 
mean that the Hawking temperature were zero for 
the models under study; 

• it is also clear that we still have the deceleration- 
acceleration transition, as we show in the plot of the 
relation for q{z) and also for q{a) in Fig. where 
our models are compared with a standard ACDM 
resulting, as said above, barely distinguishable. 

The models we have studied here involve a mixture of 
matter and the dark energy fluid which is typically the 
energy of vacuum with small modifications due to the 
variability of c and G. This means that the discussion 
of the structure formation problem (perturbation equa¬ 
tions, the formation of the structures, linear growth rate) 
is similar to those of dynamical vacuum models given in 
Ref. [43| with 7 parameter here being analogous to v pa¬ 
rameter of that reference. In fact, the models the mod¬ 
els III and IV of Ref. [43] are indistinguishable from 
ACDM while the models / and II exhibit some differ¬ 
ence what can be seen from Fig. 1 of [13], where the 
density contrast and the linear growth rate of clustering 
are shown. 


IX. CONCLUSIONS 

In this paper we extended the entropic cosmology onto 
the framework of the theories with varying gravitational 
constant G and varying speed of light c. We discussed 
the consequences of such variability onto the entropic 
force terms and the boundary terms using three differ¬ 
ent approaches which possibly relate thermodynamics, 
cosmological horizons and gravity. We started with a 
general set of the field equations which described varying 
constants entropic cosmology with a general form of the 
entropic terms. In the first approach we derived the con¬ 
tinuity equation from the first law of thermodynamics, 
Bekenstein entropy as well as Hawking temperature to 
fit the general entropic terms to this continuity equation. 
We found appropriate single-fluid accelerating cosmolog¬ 
ical solutions to these field equations. We also discussed 
the constraints on the models which come from the sec¬ 
ond law of thermodynamics. In the second approach we 
derived the entropic force for varying constants, defined 
the entropic pressure, and finally modified the continu¬ 
ity and the acceleration equations. Then, we determined 
the Friedmann equation and gave single-fluid accelerat¬ 
ing cosmological solutions as well. Finally, in the third 
approach we got gravitational Einstein field equations us¬ 
ing the heat flow through the horizon to which Bekenstein 
entropy and Hawking temperature were assigned. 

We have also examined some of the many-fluid 
(first accelerating and then decelerating) entropic mod¬ 
els against observational data (supernovae, BAO, and 
CMB). We have used data from JLA compilation of 
SDSS-H and SNLS collboration (supernovae), WiggleZ 
Dark Energy Survey and SDSS-HI Baryon Oscillation 
Spectroscopic Survey (BOSS) as well as Planck data 
(CMB). We found that the observational bound on the 
Hawking temperature coefficient 7 was much smaller 
(10“^ — 10“^) than it is usually assumed on the theo¬ 
retical basis to be of order of unity 0(1). We have also 
found that in our entropic models G should be dimin¬ 
ishing while c should be increasing with the evolution 
of the universe. Our bound on the variation of c being 
Ac/c ^ 10“® is at least one order of magnitude weaker 
than observational bound obtained from analysis of the 
quasar spectra. 
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TABLE I. Observational parameters of the entropic models under study. 
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FIG. 2. (Left.) Varying G scenario: 68 % and 95% confidence levels for q and 7 . (Right.) Varying c scenario: 68 % and 95% 
confidence levels for n and 7 . 











